:arXiv:astro-ph/9801285vl 28 Jan 1998: 



GENERAL STATISTICAL PROPERTIES OF THE 
CMB POLARIZATION FIELD 



P.D. Naselsky 

Rostov State University, Zorge 5, Rostov-Don, 344104 Russia, 
Theoretical Astrophysics Center, Juliane Maries Vej 30, 2100 Copenhagen, Denmark. 

and 
D.I.Novikov 

University of Kansas, Dept. of Physics and Astronomy, Lawrence, Kansas, 66045 
University Observatory, Juliane Maries Vej 30, 2100 Copenhagen, Denmark. 

Abstract 

The distribution of the polarization of the Cosmic Microwave Background 
(CMB) in the sky is determined by the hypothesis of random Gaussian dis- 
tribution of the primordial density perturbations. This hypotheses is well 
motivated by the inflationary cosmology. Therefore, the test of consistency 
of the statistical properties of the CMB polarization field with the Gaussian- 
ity of primordial density fluctuations is a realistic way to study the nature 
of primordial inhomogeneities in the Universe. 

This paper contains the theoretical predictions of the general statistical 
properties of the CMB polarization field. All results obtained under assump- 
tion of the Gaussian nature of the signal. We pay the special attention to the 
following two problems. First, the classification and statistics of the singular 
points of the polarization field where polarization is equal to zero. Second, 
the topology of contours of the value of the degree of polarization. We have 
investigated the percolation properties for the zones of "strong" and "weak" 
polarization. We also have calculated Minkowski functionals for the CMB 
polarization field. All results are analytical. 

Subject headings: cosmic microwave background, cosmology statistics, 
observations. 

1 Introduction 

Observations of the anisotropy and polarization of the Cosmic Microwave Back- 
ground (CMB) provide a unique information about the primordial inhomogeneity 
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of the Universe. Since detection by COBE (Smoot et al, 1992, Bennett et al, 1996) 
of the CMB anisotropy several groups have reported on high angular resolution 
observational data at the angular scales 6 ~ 1° in the vicinity of the so called 
Doppler peak in the AT/T power spectrum (Hancock, et al, 1994; Gundersen, 
1993; De Bernarais, et al., 1994; Masi, et al, 1996; Tanaka, et al., 1995; Cheng, et 
al., 1994; Netterfield et al., 1996; Scott, et al, 1996). Determination of the spec- 
trum of the primordial anisotropy on scales 6 ~ 1° will yield valuable clues to 
the formation of the large scale structure of the Universe and the most important 
parameters of the Universe: the total and baryonic densities at the present time 
(Q and f2&); the Hubble constant, ionization history etc. However, the interpre- 
tation of these experimental results as well as the comparison with the expected 
power spectra of anisotropy in different cosmological models is complicated. Fu- 
ture experiments (MAP and Planck) will construct the map of the CMB with 
high level of resolution and sensitivity. 

The temperature distribution in the sky also contains the information beyond 
power spectrum. Gaussianity of the primordial signal is well motivated by the 
inflationary cosmology and has been adopted by many authors (see for review 
Starobinsky 1982, Bardeen et al. 1983, Bardeen et al. 1986). In this case the dis- 
tribution of the CMB temperature in the sky is in the form of the two-dimensional 
scalar random Gaussian field. This field can be completely characterized by its 
power spectrum. Many authors proposed to calculate various statistical char- 
acteristics of the CMB anisotropy as tests of Gaussianity. All these techniques 
provide information beyond power spectrum: 

1. Statistic of peaks in a random Gaussian fields. Following classical pa- 
pers of Doroshkevich 1970 and Bardeen, J. M., Bond, J. R., Kaiser, N.,& Szalay, 1986 
- (BBKS)), this approach has been developed by (Bond and Efstathiou 1987, 
P.Coles 1988) for CMB anisotropy. 

2. Higher-order correlations - 3, 4 ets.(Luo and Schramm 1994, Smoot et al. 
1994, Kogut et al. 1996); 

3. Minkowski functionals as a morphological descriptors of the CMB anisotropy 
maps (Shmalzing and Gorski 1997, Winitzki and Kosowsky 1997). As it was men- 
tioned by Shmalzing and Gorski 1997, Minkowski functionals are very sensitive 
to non-Gaussianity. This approach is very effective because these functionals are 
additive with respect to the isolated regions in the sky and they have simple 
analytical form in the case of the Gaussian field. This approach can be used to 
test predictions of Gaussianity. 

4. Percolation and cluster analysis. This statistical method based on a very 
attractive idea: if the experimental signal is a sum of the primordial signal and 
non-Gaussian noise (for example foreground sources, dust emission and so on), 
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then statistical properties of the pure Gaussian signal could be distorted. This 
effect provides a basis for the investigation of the characteristics of the non- 
Gaussian noise in the experimental data. Note, that percolation has become a 
popular term among cosmologists. The percolation technique has been success- 
fully applied for investigation of the evolution of the spatial density distribution 
in the Universe due to gravitational instability (see for review Zeldovich 1982, 
Shandarin 1983, Dominik, & Shandarin 1992). For CMB anisotropy this tech- 
nique has been developed by (Naselsky and Novikov D. 1995, Novikov D. and 
Jorgensen 1996). 

Therefore, the data analysis of the CMB anisotropy can be divided into two 
parts: power spectrum estimation with subsequent cosmological parameters ex- 
traction, and investigation of the nature of the observed signal. 

There is another important characteristic of the distribution of the CMB on 
the sky: this is the CMB polarization. 

The idea that polarization provides important information about the primor- 
dial cosmic plasma was pointed out by Rees (1968). The properties of the power 
spectrum of the CMB polarization field were analyzed in for example (Basco and 
Polnarev 1979, Polnarev 1985, Bond and Efstathiou 1987, Coulson et al. 1994, 
Crittenden et al. 1995, Zaldarriaga and Harari 1995, Ng K.L. and Ng K.W., 
1995, Kosowsky 1996; Kamionkowski et al. 1996, Jungman et al. 1996, Naselsky 
and Polnarev 1987, Ng K.L. and Ng K.W. 1996, Hu and White, 1997). 

The polarization field also contains information beyond power spectrum which 
also can be used for investigation of the nature of the primordial inhomogeneity 
in the Universe. Statistical properties of the polarization field caused by Gaussian 
fluctuations was partly discussed by Bond and Efstathiou 1987, Arbuzov et al. 
1997a, Arbuzov et al. 1997b. It is important to note, that polarization contains 
more information about nature of the primordial signal than the anisotropy (the 
polarization field is a combination of two random independent Gaussian fields 
(Bond and Efstathiou 1987, while anisotropy of the CMB is only one). 

In this paper we focus attention on the general statistical properties of the 
CMB polarization field. This is not a scalar field (unlike the anisotropy) and 
can be completely described in terms of Stokes parameters - Q, U and V. Since 
Thomson scattering does not produce circular polarization, V = we can consider 

the level of polarization, which depends only on two parameters - II = — — j , 

where I is the total intensity, and P = II x I = \J Q 2 + U 2 is the polarized 
intensity. Therefore, polarization field can be described in terms of the angle 
of polarization 2ip =arctg^ and polarized intensity - P. Since polarization of 
the radiation does not have any direction (it has only the orientation (p and 
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intensity P), it cannot be formally interpreted as a vector field. Nevertheless, 
below we use the term "vector of polarization" P (so that P = \P\) for simplicity, 
taking into account that this "vector" is not directed. We assume that Q and U 
components of the auto correlated pseudo-vector P are statistically independent 
(Bond and Efstathiou 1987) and have a Gaussian distribution on the sky We 
are interested in the general statistical properties of the P distribution such as 
surface density and classification of the non-polarized points P = in the sky, 
Minkowski functionals for the value of P and the percolation of the relatively 
strongly polarized spots. 



2 Pattern of the polarization fluctuations 

In this section we discuss very specific features of the polarization pattern of the 
cosmic microwave background. All results were obtained under the assumption 
that the polarization field is the result of a random Gaussian process. We describe 
the statistical properties of the two dimensional vector field of the polarization 
P such as the surface density of the singular points P = (section 2.1), genus 
curve for the two dimensional scalar field |P| and the level of percolation through 
the relatively strongly polarized spots (section 2.2). In this section we consider 
small angular parts of the sky without loss of generality. Thus the geometry is 
approximately flat and the vector of polarization can be described in the following 
form: 

P = P x x + P y y, (1) 

where x and y are the unit vectors of the Cartesian coordinate system on the 
small angular part of the unit sphere and components P x and P y can be expressed 
in terms of Stokes parameters Q and U : 



Pcos((p) 

(2) 



P y = Psin{ip) 



Q = Pcos(2ip) 
U = Psin(2(p) 

Where ip is the orientation of polarization and: 

P x = cos(tp)Q - sin(ip)U , . 

P y = sin(ip)Q + cos(ip)U 

Therefore, components P x and P y are also independent random two-dimensional 
Gaussian fields with the same parameters as Q and U. It means, that the sta- 
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tistical properties of the vector Qx + Uy are equivalent to the properties of the 
vector P. It allows us to use the usual terms Q and U instead of P x and P y . 

2.1 Singular points of the polarization vector field. 

First, we are interested in the statistics of the singular points (xo,?/o) of the 
vector P: P(xo,2/o) = 0. This condition means, that both components Q and U 
are equal to zero in such points simultaneously: 

Q(xq, 2/o ) = U {xo,Vo) = 0. (4) 

The surface density of these points can easily be computed analytically. Points 
(xo, yo ) are the points of the intersection of the lines of zero level of the Q and U 
surfaces. The angular density of such points can be found by using the properties 
of the joint probability function for distribution of Q, U, Q\, U\, Q2, Ui- Here 
Qi and Ui are the first derivatives of Q and U respectively in the point (xo,yo): 

Of 



h{xo,yo) 

x ,y °v 



(5) 



These 6 different values are independent (Bardeen et al. 1986) for an arbitrary 
point (x, y) of the map and have zero average and the following variances: 

(Q 2 ) = (U 2 ) = al (Qf) = (Uf) = a\l% i = 1, 2. (6) 

where a Q and a\ are the spectral parameters, as they were defined by Bond and 
Efstaphiou 1987. The joint probability for these values is: 

X(Q, U, Qi, Ui)dQdUdQidUi = ■= 4 . ^ A dQdUdQidUi 

v /(27r) 6 (To<Tf 

In the vicinity of the singular point (xo, yo), the value of Q and U can be described 
by the following expression: 





(8) 

The substitution of dQdU = \ det(Qj, Ui)\dxdy and integration over dQidUi gives 
us the number density of the singular points: 

Nnp = (27t) 3 / e_( ^ +U * ) l9i' u 2 - u 2 qi\dqidui (9) 
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where q { = j£, m = ^, i = 1, 2. 

The analytical calculation of the surfase density of the singular points can be 
found in Appendix A. Here we present the main results and conclusions only, 
a. Classification of singular points 

We investigate the polarization vector field around the singular points in the 
following way. Let us imagine that a point in the vicinity of the singular point is 
moving along the lines of the vector field. In this case the investigation is similar 
to that for the singular points of linear differential equations. Following Eq.(5) 
the field in the small vicinity of the point [xq, yo), where q(xo, yo) = u(xq, yo) = 
can be described in terms of the matrix M of first derivatives of the field q and 
u, and we can consider the following equation 




where x = xq + Ax, y = yo + Ay, and 

M=( qi q2 
\ ui u 2 

The characteristic equation for Eq.(lO) is: 

X 2 -(q 1 +u 2 )X+(q 1 u 2 -q2Ui) = 0. (11) 

This equation has two roots: Ai and A2, ordered by Re(Xi) > Re(X2), and the 
classification of the singular point (xo,yo) depends on their values. 

1. If Im(Xi) = —Im(X2) / 0, then this is a focus and the vector field will spiral 
toward the point (xo,yo) (Fig.l left). 

2. If Ai and A2 are real, then the matrix M has two eigenvectors which correspond 
to different values Ai and A2, and we can consider two different cases: 

a) . A2 > - both values are positive or Ai < - both values are 
negative. This means, that the point (xq, yo) is a knot and the lines of the vector 
field tend to be aligned to the direction of the eigenvector with maximal value of 
\Xi\, i = l,2 (Fig.l middle). 

b) . Ai > and A2 < - values with opposite signs. In this case the 
point (xo,yo) is a saddle (Fig.l right). 

Singular points of different types determine the behavior of the vector field in 
their vicinities. The distribution of the singular points on the map of cosmic mi- 
crowave background polarization determines the topology of the relatively small 
polarized zones Fig.2. 
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b. Surface density of the singular points. 

Detailed calculation of the surfase density of the singular points is in the 
Appendix A. The surface density of singular points of different kinds are: 

N = 1 

N k = j§(V2-l)^ (12) 

where Nf, N k , N s are the number densities of focuses, knots and saddles re- 
spectively and r c = ^ is the correlation radius. The total number density of 
non-polarized points is: 

N np = N f + N k + N s = -L (13) 

Density of the singular points depends essentially on correlation radius (Eq. 12, 
13) and therefore on the spectral parameters o"o, o\. These spectral parameters 
depend on the spectra of polarization and on the device resolution (Bond and 
Efstathiou, 1987) (see also Fig.3). The ratios 

N f /N k = V2 + l and N f /N s = ^- (14) 

are the spectral independent constants determined only by the Gaussian nature 
of the primordial inhomogeneity in the Universe. These ratios are a characteristic 
feature of the CMB polarization vector field in the inflationary cosmology. Note, 
that for example in the two-dimensional potential vector field (V = \/a) the 
number density of foci is equal to zero since this field does not have a rotational 
component. 

The joint probability for the distribution of the eigenvalues of Ai and A2 in 
the singular points which can be either knots or saddles (not foci) is: 

F(Ai, A 2 )dAidA 2 = — ^-p-^|AiA 2 |(Ai - \ 2 )e-^-^d\ 1 d\ 2 . (15) 
4 — V2 v 71 " 

Note, that this distribution is universal for all kinds of spectra of polarization as 
well as ratios (15). 



2.2 Percolation pattern for polarization 



IPI 

Here we present our results for the Genus statistics of the value p = — . We 
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can divide the map of polarization of the CMB into two parts: regions with 
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relatively strong polarization p > po ("strongly polarized zones") and regions 
with relatively weak polarization p < po ("weakly polarized zones"). Below 
we find the value po where percolation through the "strongly polarized zones" 
changes to percolation through the "weakly polarized zones" (Fig. 3,4). Let 
us suppose that we can measure only a signal with polarized intensity p > p t , 
where pt is the threshold which determines by the sensitivity of the device. If 
we can measure only "strongly polarized" signal - p t > Po, then we can see only 
the separated polarized spots which do not percolate. Therefore, the percolation 
trough polarized zones can be reached only by the device with the sensitivity 
Pt < Po- 

The value po can be found analytically in the following way. We consider the 
value \P\ as a two-dimensional random scalar field with a Rayleigh distribution 
(Coles and Barrow, 1987). This field can be imagined as a two-dimensional 
surface in a three-dimensional space. This surface has extreme points such as 
maxima, minima, saddle points and singular points. The last ones have been 
considered in the previous subsection. The densities of maxima, minima and 
saddle points have some distributions with p: 

oo 

N m ax(p) = J n max (p')dp', 

V 
oo 

N m in(p) = J n min (p')dp' , (16) 

V 
oo 

N S ad(p) = J n sad {p')dp', 

V 

where n max (p), n m i n (p), n sa d{p) - are the number densities of maxima, min- 
ima and saddle points respectively on some interval - (p,p+dp), and N max (p), 
Nmin{p)-, N sa d{p) are the number densities of maxima, minima and saddle points 
respectively above some level p. (We note, that here saddle points are the saddle 
points of the two-dimensional surfase of p(x,y). These points are not the same 
as saddle (kind of singular points) in the previous section.) 
The definition of the Genus is: 

g(p) (p) + n min {p) - n sad (p). (17) 

The integrated Genus is then: 

oo 

G{p) = N max {p) + N mm (p) - N sad (p) = J g{p')dp'. (18) 

v 
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The level of percolation po has to be found from the condition G(po) = 0. We 
recognize that this condition does not automatically mean that po is the level of 
percolation for an arbitrary scalar field. It is well-known that for the Gaussian 
random field the percolation level corresponds to the level where Genus curve 
intersects the zero. We have checked this condition for the Rayleigh distribution 
by simulating a large number of realizations for a two-dimensional field. In the 
case of Rayleigh distribution this condition also mean that level po corresponds 
to the percolation contour. 

The detailed calculation of the Genus can be found in Appendix B. Formally 
the steps of its calculation are as follows: 

1. The value p is a combination of the independent random values q and u. The 
first and second derivatives of them are: Ui, qij, Uij, (qij = Qij/cr2,Uij = 
Uij/(T2), i = 1,2, where oi is also the spectral parameter as it was defined by 
Bond and Efstaphiou 1987: a\ = = (U^). These values obey the following 
conditions: 

p 2 = q 2 + v? 
Pi = qq-i + uui, 

IPiPj + Pij = 7(<Mj + UiUj) + qqij + uuij, 
(qu) = {q,jUj) = (qijUki) = (qu^ = (qiu) = 0, 

(qqij) = (uuij) = -%5ij, (• > 

{qiqj) = (uiUj) = \Sij, 

(qijUki) = USikSji + SuSjk + SijSki), 

^ = — 1— 

2. The joint probability F of the Gaussian distribution for the values q, qi, q^, 

Fdqdudqiduidqijduij = -j=^==e~ '4 dqdudqiduidqijduij , 
A = v x M _1 x v T 

where M is the covariance matrix and A is the quadratic form of the 12-dimensional 
vector v = (q, q i: qij,u, Ui,Uij). 

3. The substitution of p, pi, p^ in Eq.(20) from Eq.(19) and integration over 6 
variables gives us the joint probability fdpdpidpij for values p, Pi, p^ to be in the 
range from p, p i: p^ to p + dp, pi + dpi, p^ + dpij. 

4. The differential density of the extreme points obeys the equation: 

nextip) = -§ / \det(p ij )\fS(pi)S(p 2 )dp ij , (21) 



10 



where n ex t(p) is the density of the extreme points. These extreme points can be 
maxima, minima or saddle points depending on the limits of the integration over 
diiij. These limits determine by the values of tr(pij) and det{pij) of the second 
derivatives matrix (pij) (see Appendix B). 
5. The Genus curve obeys the equation: 

a 2 f 

g(p) = n ma x(p) + n min (p) - n sad (p) = -§ / det{p ij )f{p,p i = 0,pij)dpij. (22) 

<J\ J 

After integrating this, we have 

g(p) = -^(-) 2 p(p 2 -3)e- 1 r. (23) 



Air \<JoJ 
The integrated Genus curve is 

G(p) = ^y-l)e- P 4, (24) 

Condition G{p) = gives us the value of po~. 

Po = 1- (25) 
Taking into account that random value p has distribution pe 2 dp we can obtain 

that percolation through the "strongly polarized" zone when a part e 2 of the 
map is detected as a "strongly polarized". This corresponds to « 61% of the 
map's area. 

When po = in Eq.(25) we have 

»(«>) = ~ ■ (26) 

This value exactly coincides with N np in Eq. (13) with the opposite sign. These 
null-points are the non-smooth minima of the surface p. The non-smooth minima 
have not been taken into account in Eqs.(16)-(25). Therefore, the total number 
of minima per unit area is N m i nto tai = N np + Nmin(0), where N np are minima, 
if p=0 and iV m j„(0) are minima, if p^,0. Therefore the total number of extreme 
points per unit area are: 

N m axtotal -^maa;(0)) 

N sad total = N sad (0), (27) 

N m intotai = N np + Nmin(0). 
Taking into account equations (13,18,25,27) we obtain: 

Nmaxtotal ~\~ N m i n i f a i N sa( itotal 0) (^8) 

as it should be. 
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3 Minkowski functionals for CMB polarization field 



As it was mentioned above, the CMB polarization at any point of the map can be 
characterized by the orientation angle and polarized intensity - p. This intensity 
has the random Rayleigh distribution on the sky - p = \J q 2 + u 2 . Therefore, 
the value of p can be considered as a two-dimensional random Raleigh field. It 
is well-known, that two-dimensional field has only three Minkowski functionals 
which satisfy additivity and translational invariance (Minkowski 1903, Hadwiger 
1959). 

Geometrical interpretation of the Minkowsky functionals on the two-dimensional 
map is essentially easy. Analogously to the previous section, we consider polar- 
ized intensity as a two-dimensional surface in a three-dimensional space. If we 
cut this surface at the different levels pt, then the area of the map will be divided 
into two parts: the area, where polarization is above the threshold pt and the 
area, where p < pt- For a two-dimensional distribution, Minkowski functionals 
correspond to the following values: 

1. A - fraction of the area of the map, where p > pt\ 

2. L - length of the boundary between fractions, where p > pt and p < pt per 
unit area; 

3. G = N max + N m i n — N sa d - Euler characteristic (equivalent to the genus) per 
unit area. 

Therefore, threshold is the independent variable on which these functionals 
depend. The third functional has already been considered in the previous section. 

_El 

The obvious first one is e 2 . The second one can be obtained in the same way 
as it was done for the the Gaussian field (Adler 1981). Here we present the result 
without derivation: 

1 d. 

L = -p t e^- (29) 

r c 

The comparison of the Minkowski functionals for the CMB polarization field 
with these functionals for CMB anisotropy is in the Fig. 5. Functionals for 
Rayleigh distribution are equal to zero for p t < 0. The third functional should 
be described together with the number of non-polarized points (see previous 
section). These functionals can be used as a morphological descriptor of the 
CMB polarization field in a similar manner as for the CMB anizotropy (Winitski 
and Kosowski 1997). 
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4 Discussions 



In this paper we have presented calculations of the statistical properties CMB 
polarization maps. 

We believe that these statistical properties can be useful for checking the po- 
larization patterns for presence of the non-Gaussian noise (for example, confusion 
signal from sources which can have the same spectral parameters as the polariza- 
tion of the CMB). If an observational signal is free from non-Gaussian noise and 
is Gaussian itself (due to inflation) then the topological approach is not neces- 
sary, because the correlation function or equivalent its power spectrum contains 
all information about the polarization signal. On the other hand, if the signal is 
a sum of polarization of the CMB (which is Gaussian) and unresolved foreground 
sources (which are non-Gaussian), then a detailed topological picture of the po- 
larization field around the non-polarization points will be distorted comparative 
the predictions of the theory for the Gaussian distributions. 

On the other hand, the investigation of the nature of the primordial polar- 
ized signal is a test of the inflationary model of the evolution of the Universe. 
Therefore, the investigation of the Minkowski functionals together with the non- 
polarized points on the observational data and comparison with the theoretical 
predictions can be used as the test on Gaussianity of the primordial inhomogene- 
ity 

We would like to emphasize that the regions with strong polarization will be 
detected easier than the regions with weak polarization. As we demonstrated 
in the paper (see section 2.2) these region occupy an essential part of a whole 
map. From this point of view it is interesting to study the statistical properties 
of these regions. It is worth also mentioning that it is interesting to investigate 
the dependence of the spectral parameters of polarization in various cosmological 
models on the resolution of the detector and related statistical properties of the 
maps of polarization of the CMB. It is also very interesting to stady the cross- 
correlations between anisotropy and polarization on the sky map and make some 
theoretical predictions from the geometrical point of view. These quations will 
be considered in a separate paper. 
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Appendix A 
Surface density of the singular points 

Below we describe the density of foci, knots and saddles. The total density of 
singular (non-polarized) points is: 



4 f a\\ f -4. 
N np = -j—^ y—j J \q\u 2 - uiq 2 \e 2 dqidu h (Al) 



where 

A= J2 2g? + 2«?, z = 1,2. 

i=l,2 

The substitution 

q\ = tj(x + Rcos ip), u\ = \ (w + R sin ip) , 
(?2 = \ (R sin ip — w) , U2 = \ (x — R cos 99) , 

and integration over d(p gives us: 

1 1 /crA 2 



(A2) 



j e -^+B?+w^ x 2 +w 2_ R 2\ R dRdxdw (A3) 



The next substitution b = R 2 — w 2 allows us to rewrite (A3) in the following 
form: 



Nnp= OV e ' hx2dx I e ~ w2dw I \ x2 ~ b \ e ~ idb ^) 



where — w 2 < b < 00, —00 < x < +00, —00 < w < +00. In terms of x and b the 
eigenvalues of the matrix 

M = { qi q2 \ {Ah) 

Ul «2 
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are 



\ 1 , 2 = ±{x±Vb). (A6) 

From (A4, A6) we can obtain the density of focuses, saddles and knots: 
-w 2 < b < - foci 

< b < x 2 - knots (A7) 
x 2 < b < oo — saddles 

Using (A4), (A7) we obtain 

16-7T ^ (To ) 

^ = w(2- 2 ^2)fe) 2 (AS) 

Using (A4, A6-A8) the joint probability for values Ai, A2 in the peculiar points 
which can be knots or saddles (not focuses) is 

P(Ai, X 2 )dX l d\ 2 = 8 ——\\ l \ 2 \(\ l - \ 2 )e- x '- x 2dX 1 d\ 2 , X 1 > A 2 . (A9) 

(4 - V^JVTT 



Appendix B 

Genus curve 

In this appendix we obtain the differential and integrated Genus curve for the 
two-dimensional random Rayleigh field. 

According to section 2.2 the value P = \/q 2 + u 2 is the non-linear combination 
of two different independent random Gaussian fields q and u. Equation (17) for 
the joint probability distribution of the values q, u, qi, Ui qij, Uij, i,j = 1,2 
contains the quadratic form A and det M, where M is a covariance matrix: 

A = q 2 +u 2 + 2(q 2 + q\ + u\ + u 2 ) + ^+^+^f+^+^+^) 2 + 

+2(9ii - q 2 2? + 2(«n - u 22 ) 2 + 8q 2 2 + 8u 2 2 , {Bl) 

detM = 2- 12 (l- 7 2 ) 2 . 
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The substitutions in Eq.(17) 

= psin (p, 

— ri ■ cm in — 11 ■ pnc in 

(B2) 



q = pcosip, u = psm(p, 

Pi = qi cos if + Ui sin if, hi = q, L sin ip — m cos ip, 
lij = qij cos ip + sin ip, lij = qij sin <p — Uij cos 99, 



and integration over dkdkj gives us the joint probability for the distribution of 
the values p, pj, /«, Z^: 

32 A 

X(p,pi,li, lij)dpdpidlidlij = (2^)7/2^ _ ^ 2 y /2 pe~~ dpdpidkdkj i,j = 1, 2 (B3) 
Following Eqs.(16) and (B2) we can get the expression: 

PPij = Pkj+lklj; {B4) 
using Eqs.(20), (B3) and (B4) we obtain 

32 f A 

9(P) = ( 27 r)7/2(i_ 7 2)i/2 J (Pnf22 " P&pe—dlidpij i,j = 1, 2 

A = P 2 + 2(Z 2 + Z 2 2 ) + ^ + ^ 7 7Q)2 +2 ( Pll -P22-7fr) 2 + 8(p 12 - 7 c) 2 (S5) 

a= l l±%-p, 6=i— 3 C = M 
p p p 

the integration over dpjj, dZj gives us differential Genus curve 

5b) = ^(^)%(p 2 -3)e^ (£6) 
The integrated curve is 

G(p) = / 9(P'W = ^ (^) (P 2 -I)e^. (£7) 
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Figure 1: Classification of singular points. Field of polarization around singular 
points p = 0. Dashed lines show the direction of the pseudo- vector P (but not 
its value). Left - focus, middle - knot, right - saddle. 
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Figure 2: Simulated map 2 x 2 of the CMB polarization field for the scale- 
invariant adiabatic CDM model with Q = 1, Qf, = 0.03, h=0.75 with smoothing 
angle 5 arcmin (FWHM). The simulation technique for small parts of the sky 
and spectrum for simulations are from Bond and Efstaphiou 1987. Left 
- polarization field. The length of eacBOvector is proportional to the degree of 
polarization and the orientation gives the plane of polarization. For visual clarity, 
we only use 50 x 50 vectors. Right - the same as left, but we plot only the 

orientation of polarization in the vicinity of non-polarized points (solid lines). 
This map containes 7 non-polarized points - 2 foci, 1 knot and 4 saddles. 
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Figure 3: Simulated map 10 x 10 of CMB polarization for the same model as 
in the fig. 2. Dashed area corresponds to the regions with polarization degree 
p > pt- Solid lines are the boundary between regions with p > p t and p < p t . 
Circles, threangles and stars are foci, knots and saddles correspondingly. This 
map containes 13 fosi, 6 knots and 19 saHdles. 




Figure 4: The same as in fig. 3, but without non-polarized points. We plot 
area, where p > pt for different values of pt- pt = 2, 1.5, 1, 0.5. Spots with p > p t 
percolate then p t = 1, which corresponds to the e -1 / 2 ~ 61% of the maps area. 
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Figure 5: Minkowski functionals for CMB polarization (solid lines) and 
anisotropy (dashed lines). Threeshold is given in the units of uq for polariza- 
tion and in the units of y/ ((AT) 2 ) for anisotropy. 
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